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1 Introduction 



It is well known that in classical mechanics the knowledge of all first integrals 
of motion of a given problem is equivalent to finding its complete solution. 
Nowadays the search for first integrals has assumed an increasing impor- 
tance in the determination of the integrability of a dynamical system. It is 
extremely important to know if a non-linear dynamical system will present 
chaotic behavior in some regions of the phase space. The notion of integrabil- 
ity is related to the existence of first integrals of motion. Several methods of 
finding first integrals are available in the literature for example, Lie's method 
KJ, Noether's theorem 0, or the direct method ||. Even if not all first in- 
tegrals of motion associated with the problem at hand are found, it may 
happen that the ones which are obtained contribute to the discovery of the 
solution we are seeking for. Nevertheless, if we do find the solution we are 
after by solving the equations of motion in a straightforward way, it still may 
be profitable to look for additional constants of motion. Such is the case of 
the Kepler problem where the knowledge of the Laplace-Runge-Lenz vector 
||, |J allows us to obtain the orbit in a simple way. 

Of the inexhaustible wealth of problems which we can find in classical 
mechanics one of the most aesthetically appealing and important is the cen- 
tral field problem. Energy and angular momentum associated with this type 
of field are well known conserved quantities. However, other vector and ten- 
sor conserved quantities have been associated with some particular central 
fields. The Laplace-Runge-Lenz vector is a vector first integral of motion for 
the Kepler problem; the Fradkin tensor J7J is conserved for the case of the 
harmonic oscillator and for any central field it is possible to find a vector first 
integral of motion as was shown in Q. In the general case these additional 
integrals of motion turn out to be complicated functions of the position r and 
linear momentum p of the particle probing the central field. When orbits 
are closed and degenerated with respect to the mechanical energy, however, 
we should expect these additional constant of motion to be simple function 
of r and p. In this article we wish to exploit further this line of reasoning by 
determining the existence of such additional vector first integrals of motion 
for the time-dependent Kepler and isotropic harmonic oscillator problems. 
In particular, we will show that for the time-dependent Kepler problem the 
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existence of a vector constant of motion coupled to a simple transformation 
of variables turns the problem easily integrable. 

The structure of this paper goes as follows: in section 2 we establish the 
conditions which guarantee the existence of a vector first integral for a general 
central force field. In section 3 we put the method to test by rederiving some 
well known results such as the conservation of angular momentum in an 
arbitrary central field, the conservation of the Laplace-Runge-Lenz vector for 
the Kepler problem, and the conservation of the Fradkin tensor fixing en route 
these specific fields to which they correspond. In section 4 we consider the 
time-dependent case establishing generalizations of the examples considered 
before and presenting new ones. In section 5 we show that the existence of 
a vector first integral enable us to find the orbits of a test particle. This 
is accomplished for the case of harmonic oscillator, and the time-dependent 
Kepler problem. Also the period of the time-dependent Kepler problem is 
obtained. Finally, section 6 is reserved for final comments. 

2 Constructing vector constants of motion 

The force f (r, t) acting on a test particle moving in a central but otherwise 
arbitrary and possible time- dependent field of force g(r,t) can be written as 



where r = r (t) is the position vector with respect to the center of force, r is 
its magnitude, and t is the time. To this test particle we assume that it is 
possible to associate a vector j which in principle can be written in the form 



where p = p (t) := mi (t) is the linear momentum, m is the reduced mass 
and A, B are arbitrary scalar functions of p, r and t. Taking the total time 
derivative of (0) and making use of Newton's second law of motion we readily 
obtain 



If we assume that j is a constant of motion it follows that the functions A 
and B must satisfy 




(1) 




(2) 




(3) 
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(4) 



dt m 

Eliminating B between (f|) and (jE)) we obtain 




(5) 



m —-gA = 0. 
It follows from (|5|) that j can be written in the form 



(6) 



. „ dA 
j = A p — rn—— r. 



(7) 



Therefore, since (||) is equivalent to both (HI) and (^) if the field g(r,t) is 
known any solution of (@) will yield a vector constant of motion of the form 
given by (|7|). Equation (]|), however, is a differential equation whose solution 
may turn out to be a hard task to accomplish. Nevertheless, we can make 
progress if instead of trying to tackle it directly we make plausible guesses 
concerning A thereby linking j to specific forms of the field g(r,t). This 
procedure is tantamount to answering the following question: given j what 
type of central field will admit it as a constant of motion? The answer is 
given in the next section. 

3 Simple examples 

With r, p and a unit constant vector u we can construct the following scalars: 
u ■ r, up and r • p. Other possibilities will be considered later on. For 
the moment let us consider some simple possibilities for the scalar function 



Consider first A(p, r, t) — u • r. It is immediately seen that this choice for 
A satisfies @ for any function g(r,t). The constant vector j reads 



and it can be related to the angular momentum L = r x p as follows. Firstly 
we recast (§) into the form 



^(P,r, 



t). 



j =(u r) p-(u p) r, 



(8) 



J = 



M ■ u 



(9) 
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where M = pr — rp. Since u is a constant vector we conclude that the 
constancy of j is equivalent to the constancy of M whose components are 
Mjk = PjXk — XjPk, with i, j, k = 1,2,3. The antisymmetrical tensor M 
is closely related to angular momentum L of the test particle. In fact, it 
can be easily shown that — 2L, = e^M^, where Lj is the i-th angular mo- 
mentum component and e^k is the usual permutation symbol or Levi-Civita 
density. Therefore, this simple choice for A leads to conservation of angular 
momentum for motion under a central arbitrary field g{r,t). 

Consider now the choice A(p, r, t) = u • p. Then making use of Newton's 
second law it follows that (|7|) is satisfied if we find a solution to 

^Mu-r = 0. (10) 
dt 

For arbitrary values of u ■ r we can find a solution to (|TI]) if and only if 
g(r, t)= 0, or g = go where go is a constant. In this case we can write 

j = (u • P) P - 9o (u • r) r. (11) 

If we choose the constant to be equal to — k then the central force field 
will correspond to an isotropic harmonic oscillator, f = —k r. As before, ( |TT| ) 
can be recasted into the form 



j = 2mF-u, (12) 

where F is given by 

F=PP + ^. (13) 

2m 2 K J 

The tensor F is symmetrical and is known as the Fradkin tensor [|7| . Finally, 
consider A = r • p. For this choice of A (0) yields 

1 dg 3 dr 

— r + — r = °> ( 14 ) 
g at r dt 

where we have made use of (|l|) and also of the fact that dr/dt and r are per- 
pendicular vectors. Equation ( 14|) can be easily integrated if g is considered 



to be a function of the radial distance r only. If this is the case we obtain 
the Kepler field g(r) = —k/r 3 . The constant vector j is then given by 

r 

j = (r ■ p) p — (p • p) r + mk-, (15) 

r 
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Making use of a well known vector identity we can recast (fHS| ) into the form, 



j = L x p — mk-. 

r 



(16) 



Therefore j can be equaled to minus the Laplace-Runge-Lenz vector A. From 
( JT5| ) and the condition j • r = —A ■ r = the allowed orbits for the Kepler 
problem can be obtained in a simple way, see for example ||. 

4 Time-dependent fields 

Let us now consider time-dependent central force fields for which we can build 
more general vector first integrals of motion. As with the time-independent 
case there are of course several possibilities when it comes to the choice of a 
function A for a time- dependent central field. Here is one 



A = (p{t) r • p + ip(t) r ■ r. 
Evaluating the second derivative of (|T7|) we obtain 

d 2 A (d 2 



dt 2 



dt* 



4# 4 dip\ 
H 1 — r ■ p 



m 



m dt 



(17) 



2^t g + 0^ + + 2 ^ 9 \ r . r 
dt dt dt 2 m J 



+ 



m- 



m z 



dt 



+ ip p • p. 



(18) 



Where we have made use of (|T]). If we impose the additional condition 



ip + m 



dt 



0. 



(19) 



we eliminate the quadratic term in p. With the condition given by ( |19| ) 
we can substitute for ip in (|17D and ( pT8|) and take the results into (§) thus 
obtaining 



3 {-^J + 



d r 

~dt 



dg_ 

' dt 



m 



d?(p d(p 



r • r 



0. 



Equation ( p0[) can be rewritten as 



3 
2 



d 2 (j) \ d(r 2 ) 



dt 2 



dt 



d 2 



d_ (_ 

dt \ m ~dl 2 



+ 



(20) 



(21) 
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and easily integrated so as to yield 



m d 2 d> C 

9(r,t) = --^ + —. (22) 



dt 2 0r 3 

The vector first integral of motion associated with is 



( , d(j> \ ( dcf) . mC\ /nnS 
j = I (pr ■ p - m ~^ v - r P + I m "^T r " P ~ ' P —J r > I 23 ) 

which can be simplified and written in the form 

, 9t d fr\ mC r .„ JS 

j =^ Lx sU)-— • (24) 

where we have made used of the fact the angular momentum is constant for 
any arbitrary central field whether it is time-independent or not. If in ( |24| ) 
we set = 1 and C^O, then from fl22|) we see that g (r, t) is the Kepler field 
and j is minus the Laplace-Runge-Lenz vector as before; the scalar function 
A (r, p, t) reduces to r ■ p which we have already employed in section 3. If we 
set (m/<j)) d 2 (j)/dt 2 = —k(t), that is, if <fr is an arbitrary function of the time, 
and also C = 0, we have the time-dependent isotropic harmonic oscillator 
field, F (r) = —k(t) r. In this case j is equal to the first term on the R.H.S. 
of (P^l) . If (m/<j))d 2 <f)/dt 2 = —k and C = 0, we have the time-independent 
isotopic harmonic oscillator field but this time j is not the same vector as the 
one we have obtained before. The reason for this is our choice ( fPTD for the 
scalar function A (r, p, t) which is not reducible to the form u • p employed 
previously. 

As a last example let us consider again the time-dependent isotropic har- 
monic oscillator and show how it is possible to generalize the Fradkin tensor 
for this case. Let the function A (r, p, t) be written as 

A = <f>(t)u-T + if)(t)u-p. (25) 
The first and the second derivative of A read 

dA d(p „ „ dib „ , A ,_„ N 

— = — u-r+— u-p + — u-p + ^u-r, (26) 
dt dt m dt 
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and 
d 2 A 



dt 2 



dt 2 dt dt m 



Taking fl27p into (Bp we obtain the condition 



(27) 



(d 2 (j) dip dg \ „ f d<f> d 2 ^ 

m U^ +2 ^ + d^ u ' r+ [ 2 n +m ^ 



up 



0. 



(28) 



Imposing the additional condition 



28T) becomes 



d<f) d 2 ^ 
2— + m—— 
dt dt 2 



dt 3 dt dt 



(29) 



(30) 



We can solve (|30"|) thoroughly if g (r, t) is a function of the time t only. In 
this CcLS6 ; clS before, we end up with the time-dependent isotropic harmonic 
oscillator. The vector j associated with (|25|) can be obtained as follows: first 
we integrate ( f29"D thus obtaining 



mdtp 

where C is an integration constant. Then making use of 
we arrive at 



(31) 

(0) and (EH) 



m d^ 



+ 



2 

'm 2 d 2 ip \ „ fmdip \ ^ 

u-r- + C u-p 



(32) 



or in terms of components 



Ji -Mi ? W 



(33) 



where .FL is defined by 
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( mdtb \ 

F ij = l-^-^ + CjPiXj+ippiPj (34) 

/ m 2 d 2 ip \ fmdijj \ 

+ - mgiP ) x iXj - + Cj x lPj . 

The constant C in ( [34]) can be made zero without loss of generality. A 
generalized Fradkin tensor can now be defined by 



ij = l~ ~2~dl ) PiXj ^ PiPj \T1F ~ mg ^ ) XiXj ~ ~2~dl XiVj ' ^' 

From (^) we can read out the diagonal components of the generalized Frad- 
kin tensor 

2 d ip dxi 2 ( m2 d 2 ^ \ 2 
F, = ~m— Xl — +1 p Pi + ^ T — - mg j Xi . (36) 

It is not hard to see that the trace of this generalized Fradkin tensor F (r, p, t) 
becomes the energy of the particle when g (r, t) is a constant field. 



5 Obtaining explicit solutions: An alterna- 
tive way 

Now we wish to show how to take advantage of the vector constant j to obtain 
the solution for the Kepler and the isotropic harmonic oscillator potentials. 
But firstly we must establish some very general relationships between the 
sought for solution r (t) and A (r, p, t) and j. First notice that (|T^) can be 
recasted into the form 

(37) 

where it must be kept in mind that A (r, p, t) satisfies (^|). As we have seen 
before in specific examples the form of the vector constant j depends on the 
force acting on the particle. Integrating ([37]) we readily obtain 



2 d 



j = m A(r,p,t) — 



A(r,p,t) 
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A(r,p,t) A(r,p,0) m Jo A(r,p,r) 

Equation (|38|) can be given a simple but interesting geometrical interpreta- 
tion. Assume that the initial conditions r (0) and p (0) are known and there- 
fore the function A (r, p, 0) can be determined. The vector r (0) j A (r, p, 0)is 
therefore a constant and completely determined vector. As time increases, 



the R.H.S. of fl38|) increases. The vector on the right side of (|38|) though 
varying in time has a fixed direction which is determined by j. Therefore, 
r it) /A(r,p,t) must increase in order to close the triangle whose sides are 
the three vectors involved in fl38l). If the orbit is unlimited then it is easy 
to see that the following property ensues: there is an asymptote if in the 
limit t — > oo the definite integral Jq is constant. On the other hand, 
if the orbit is limited, but not necessarily closed, there will be a position 
vector r whose direction is parallel to that of the vector j at the instant 
t*. If the length of the position vector r is finite, we can conclude that 
at the same instant t* the function A (r, p, 0) must be zero. Thus, we can 
see that the vector r (t) /A(r,p,t) must be reversed at this instant and its 
evolution is determined by the fact that its end is on the straight line that 
contains j. For t = t* + e, where e is a positive infinitesimal number, the 
vector r (i* + e) /A (r, p, t* + e) changes its direction abruptly, so to speak, 
as shown in the figure, hence in the transition A (r, p, £*) — > A (r, p, t* + e) 
the scalar function must change its sign. 

Let us obtain the solution r (t) for the case of the isotropic harmonic 
oscillator. A particular solution of ( |T7D for g — —k, where k is the elastic 
constant is given by 

A(t) = cos{ut) , (39) 



where u = y— is the angular frequency. The solution given by (|39|) allows 
us to write 

r( 0) = -EM/"^l^ (40) 



cos(o;£) m Jo cos 2 (lot) 

The integral can be readily performed and after some simplifications we fi- 
nally obtain 

r (t) = cos uot r (0) + — sin ut p (0) . (41) 



10 



Therefore we have obtained the solution of the time-independent harmonic 
oscillator in an alternative way from the knowledge of the initial conditions 
as it should be. Notice that the general solution A(t) = A (0) cos (out + 9) 
would lead to the same general result. Another possible solution in the case 
of the time-independent isotropic harmonic oscillator is given by 

A(r,p,t) = u.p(t), (42) 

as can be shown by substituting this solution into (|6]). This solution shows 
that the trajectories have no asymptotes. 

Let us now show how we can obtain the orbits in the case of the time- 
dependent Kepler problem. Let us begin by rewriting ([J7J) in polar coordi- 
nates on the plane. In these coordinates the angular momentum conservation 
law is written in the form 

I = mr 2 — (43) 
dt v ; 

and this allows to rewrite (071) as 



A A 

Introducing the unitary vectors r and 9 we can write the above equation as 



d fA\ a A a 

r +- e 



d9 \ r J r 

A A 

The components of the vector j in the direction of r and 9 are given by 



(45) 



, a A j . 

j- 9= I- , (46) 
r 



and 



a 



d (A\ 



j ' r = - l To \-) ' (47) 

Equation ([IT]) can be obtained from ( pf) and therefore it is redundant. In 
section 4 we determined a generalized Laplace-Runge-Lenz vector for the 
time-dependent Kepler problem. The scalar function A (r, p, t) associated 
with this vector was found to be 
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A = <p[t) r • p — m—r 

(JjL 



(48) 



Making use of (13) we can rewrite the linear momentum as a function of 8 
as follows 



dr d6 I dr 
p = m dt dt = ^~d~6 



Taking (|49| ) into ( ^8|) and considering A as a function of 9 we obtain 



r d9 \r 



(49) 



(50) 



Equations. (|46|) and (|^) lead to 



dd [r 



y± = j_ 
i 2 p 



sin 



a 



(51) 



where a is the angle between j and the OX axis (see figure 2). In order 
to integrate (|5l|) we assume that the initial conditions at t = are known 



vector functions, i.e. 



and 



r(0) = r ; 



(52) 



P(0) = Po • 

In terms of polar coordinates these initial conditions are written as 



(53) 



r(9 ) = r 



and making use of (49) 



dr^ 
d9 



= 00 



Upon integrating ( pTj ) we find 



t = ll _|_ 2 
r r I 



+ — [cos(6 l - at) — cos( 



-a)] 



(54) 



(55) 



(56) 
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For the usual time- independent Kepler problem, 0=1 and in this case (SB) 
takes the form 



I j 

+ — [cos(^ — a) — cos(# — a)] 



r 



(57) 



The scalar product between j as given by (16) and ro permit us to eliminate 
008(6*0 — a) and leads to the usual orbit equation 



mC 
1?~ 



mC 



cos(# — a) 



(58) 



If we define a new position vector r' according to 



and redefine our time parameter according to 

dt : = T2 ' 



(59) 



(60) 



we can recast the equation of motion for the time-dependent Kepler problem, 
namely 



rn 



'md 2 (h C 
+ 



(j) dt 2 (f)r 3 



(61) 



into a simpler form. According to (^) and (^) the velocity and the accel- 
eration transform in the following way 



dr 

dt 



1 dr' 



and 



dt 2 



J 2 



+ 



1 d 2 r' 



(62) 



(63) 



dt 2 (j) 3 dt' 2 

where we have taken advantage of the fact that r, r' and (ft can be considered 
as functions of t or t' . Making use of (|63D the equation of motion (|B~lD can 
be written as 



m 



dt 12 



c 

VA3 1 



(64) 
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Equation (|5"3] ) corresponds to the usual time-independent Kepler problem 
whose solution is given by (|58|). Equations (|59|) and (^Cj) show that the 
open solutions of (|54]) are transformed into the open solutions of fl61"|) with 
the same angular size and that closed solutions of (|64] ) are associated with 
spiraling solutions of (^l|). The period of the orbit of (|64"D is related to the 
time interval that the spiraling particle takes to cross a fixed straight line. 
Representing this time interval by T we have 

As an application of the above remarks suppose we are looking for the form 
of the function which yields a circular orbit with radius R as a solution 
to fl61|)? To find this function we see from ( |B"T1) that we have to solve the 
following equation differential equation 

+ ZTBif^ = ZTTxv ( 66 ) 



The solution is 



dt 2 mR 3 mR 3 



0(t) = l + <f) cos{ujt + (3) , (67) 



where 0o is a constant and oj = y and (3 an arbitrary phase angle . The 
constant 0o may be chosen so that the transformed solution will be a given 
ellipse as we show below. Equation (^) leads to 

A = R- 1 [1 + O cos (tot + 13) ] . 

Comparing with we obtain 



and 



m \C\ 

The period of this circular orbit is given 



R = , (68) 

m C 



(69) 
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2tt 

T = — = 2tt, 
uj 



(70) 



and using (|65|) we get 



^ fij- dt 2n 1 . , 

T = = (71) 

7o (l + ecoscut) 2 w (l-e 2 )i 

where e = |^>o| is the eccentricity. Making use of 1 — e 2 = R/a, where a is 
the major semi-axis we finally obtain the orbital period 



T = 27 ra^--. (72) 
To conclude consider the total mechanical energy associated with (|64l) 

£ = — + — = const. (73) 
2 m r' K ' 

Since p' and p are related by 

p' = (j)(t) p -m<j)(t) r (74) 
and r' and r by ( |5T?D we easily obtain 

E = 2^- 2 ^ r ' p+ UJ T + (75) 

which is a conserved quantity and can be interpreted as a generalization of 
the energy of the particle under the action of a time-dependent Kepler field. 



6 Laplace- Runge-Lenz type of vector constants 
for arbitrary central fields 

Equation ( f22|) determines a time-dependent Kepler field g (r, t) , where the 
variables r and t are independent. If, however, we consider the orbit equa- 
tion r (t) we can eliminate the time variable and define the function g (r, t (r)) 
which can be understood as an arbitrary function of r. For the sake of sim- 
plicity we denote this function by g (r) . The function (t) which transforms 
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the Kepler problem when understood as a function of r transforms the Kepler 
field in an arbitrary central field. Let us write Eq. (p2|) in the form 



dt 2 

and consider the transformation 

d 2 6 d 2 c 



d 2 <f) . . ± C 
m^-r. — g (r, t) (p H — - 







' dr\ d<p d 2 r 
dr dt 2 



dt 2 dr 2 \dt / 

Energy conservation and the equation of motion allow us to write 

r 2 



(76) 



(77) 



and 



(dr^ 


2 


2 




) = 


m 




d 2 r 


i 




M 2 





E-V(r) 



2mr 2 



rg (r) 



ni 



+ 



m 2 r 3 



(78) 



(79) 



where E is the energy of the particle and / its angular momentum. Taking 
these three last equations into account Eq. ( |75| ) reads now 



E-V(r) 



/ 2 



2mr s 



d 2 



dr^ 



+ 



rg (r) + 



I 1 



m 2 r 3 



^-,(r)0+^ = O (80) 



Equation (|80D permit us to determine the function (j> (r) for any potential 
V(r). Thus, we conclude that a central field problem can be transformed 
into a time-dependent Kepler problem. When g (r) describes the Kepler field 
the solution of ( |30D is <p ( r ) — 1- Sometimes it is convenient to perform a 
second change of variables by defining the transformation 



= ijj 



mC 



Then Eq. 



) becomes 
I 2 



E-V(r)- 



2mr 2 



d 2 *P 
dr 2 



+ 



rg (r) 



rn 



+ 



m 2 r 3 



dip 
dr 



(81) 



-<7(r)V = (82) 



As an example consider V (r) = k/r. Then the solution of (|32D is simply 

ijj ( r ) = ci (kmr + Z 2 ) + c 2 Vl 2 + 2kmr - 2mEr 2 (83) 
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7 Conclusions 



In this paper we have outlined a simple and effective method for treating 
problems related with time-dependent and time independent central force 
fields. In particular we have dealt with the Kepler problem and the isotropic 
harmonic oscillator fields. We have been able to rederive some known results 
from an original point of view and generalize others. The central force field 
has been discussed in the literature from many points of view. The difficulty 
in finding vector constants of motion for central fields stem from the fact that 
in general orbits for these type of problem are not closed, therefore any new 
ways to attack those problems are welcome.. In our method this difficulty is 
transferred, so to speak, to the obtention for each possible central field, which 
can be time-dependent or not, of a certain scalar function of the position, 
linear momentum and time. For a given central field this scalar function is 
a solution of @. In the general case, the obtention of the scalar function is 
a difficult task. Judicious guesses, however, facilitate the search for solution 
of (^|) and this is what we have done here. 
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